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Abstract 

Recent three-dimensional radiation hydrodynamic sim- 
ulations by Wcde meyer et al. 2004| suggest that the solar 
chromosphere is highly structured in space and time on 
scales of only 1000 km and 20-25 sec, resp.. The resulting 
pattern consists of a network of hot gas and enclosed cool 
regions which are due to the propagation and interaction 
of shock fronts. In contrast to many other diagnostics, the 
radio continuum at millimeter wavelengths is formed in 
LTE, and provides a rather direct measure of the thermal 
structure. It thus facilitates the comparison between nu- 
merical model and observation. While the involved time 
and length scales are not accessible with todays equip- 
ment for that wavelength range, the next generation of 
instruments, such as the Atacama Large Millimeter Array 
(ALMA) , will provide a big step towards the required res- 
olution. Here we present results of radiative transfer calcu- 
lations at mm and sub-mm wavelengths with emphasis on 
spatial and temporal resolution which are crucial for the 
ongoing discussion about the chromospheric temperature 
structure. 
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1. Introduction 

The structure of the solar chromosphere, and in partic- 
ular of internetwork regions, is still an issue of intensive 
scientific debate. On the one hand UV diagnostics imply 
a hot layer with only small temperature fluctuations, as 
expressed in static semi-empirical 1-D models like, e.g., 
FAL A, VAL C IjFontenla et al. 1993llVernazza et al. 1981|l . 
These models feature a temperature minimum and above 
a smooth increase towards high coronal values. In contrast 
the observed amount of carbon monoxide requires cool gas 
to be formed. Appropriate models (see, e.g., |Ayres 2 002) 
thus show very low temperatures at heights where the 
aforementioned class of models predicts much higher val- 
ues. It should be noted that recent works IjAsensio Ramos et 
lUitenbroek 2000l|Wedemeyer et al. 2004| ) suggest that CO 
is mostly located in the photosphere and low chromo- 
sphere, and its presence thus does not necessarily indicates 
low temperatures in the higher layers. 



The different models are very elaborate and explain the 
observations, which they were constructed for, very well. 
However, both type of models cannot provide a complete 
description of the solar chromosphere since they fail to 
explain the entirety of all available observations at the 
same time. 

A solution might be a very inhomogeneous and dy- 
namic chromosphere which does exhibit cool and hot tem- 
peratures next to each other. Examples are the recent 3D 
radiation hydrodynamic simulation by Wedemeyer et al. 
(2004, hereafter Paper I) and the pioneering 1-D work 
done by Carlsson & Stein (e.g., 1995). Both show a very 
dynamic picture produced by propagating shock waves. 
Moreover, the mentioned 3D model reveals a network-like 
pattern of hot and cool gas on spatial scales which are of 
the same size as the underlying granulation pattern. 

In order to confirm the structure to be inhomogeneous 
and dynamic as implied by these models, detailed com- 
parisons with observations are mandatory. But still such 
direct comparisons are hampered by various difficulties. 
One major problem, e.g. in the UV, is the breakdown of 
the LTE (local thermodynamic equilibrium) assumption 
at chromospheric heights, causing the translation of tem- 
perature into intensity to be a rather involved process. 
Hence, the production of synthetic intensity images from 
numerical models, which can be compared with observa- 
tions, requires detailed multi-dimensional non-LTE radia- 
tive transfer calculations. Such calculations are hardly fea- 
sible at present, since they require large computational 
resources and efforts. 

The situation is physically simpler for the radio contin- 
uum at millimeter and sub-millimeter wavelengths which 
can be exploited as a linear measure of gas temperature, 
owing to the validity of LTE in this wavelength range. The 
necessary radiative transfer calculations can be done com- 
paratively easily. Hence, (sub-)millimeter wavelengths of- 
fer a convenient way to compare numerical models and ob- 
servations as demonstrated in a recent work by |Loukitcheva et al. ( 
Unfortunately, no appropriate observations are available 
so far since the spatial resolution of the existing instru- 
aL c 2S93~ or this wavelength domain is too low, rendering 
the desired granular scales unaccessible for now. 

The situation will substantially improve in the future, 
when the Atacama Large Millimeter Array (ALMA) will 
commence operation (~ 2011). This instrument will pro- 
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vide high spatial and temporal resolution, allowing to ob- 
serve the small-scale structure of the solar chromosphere. 

Here we present images of emergent continuum inten- 
sity at millimeter wavelengths which are calculated from 
the recent three-dimensional radiation hydrodynamics sim- 
ulation by |Wedemeyer et al. (2004) Since the simulation 
does not include magnetic fields, the presented analysis 
refers to non-magnetic internetwork regions only. Our aim 
is to predict how those areas would be seen by instruments 
like ALMA. 

A detailed paper with supplementary calculations is in 
preparation. 

2. Atacama Large Millimeter Array 

The Atacama Large Millimeter Array (ALMA) will consist 
of 64 telescopes with diameters of 12 m each, setting up 
2016 baselines. The full array will start operation in 2011 
in the Chilean Andes. There will be ten frequency bands 
which cover a total range from 31.3 GHz (A = 9.58 mm) 
to 950 GHz (A = 0.32 mm). According to |Bastian (2002)[ 
ALMA will provide an angular resolution of 0."015 to l."4, 
depending on antenna configuration. This corresponds to 
only ~ 10 km to ~ 1000 km on the Sun. 

However, ALMA will not sample both these small spa- 
tial scales and the larger scales needed to reconstruct the 
entire field of view (e.g., 21" at A = 1 mm or about the 
size of the interior of an internetwork region) simultane- 
ously. Rather, the effective resolution may be limited by 
the number of pieces of information measured instanta- 
neously: 2016 resolution elements (due to the number of 
baselines) correspond to an effective resolution of about 
0'.'5 at A = 1 mm. 

A major advantage of ALMA for solar research lies 
in the properties of its wavelength range because (i) gas 
temperature linearly translates into emergent continuum 
intensity and (ii) synthetic intensity images for compari- 
son can be calculated comparatively easily. 



3. The hydrodynamic model 

Here we use the numerical 3D model described in Paper I 
which has been calculated with the radiation hydrody- 
namics code C0 5 BOLD flFreytag et al. 2002| ). Magnetic 
fields are not included, restricting the model to internet- 
work regions. The resolution of the model atmosphere is 
40 km in horizontal (x, y) and 12 km in vertical direc- 
tion (z). The top of the model is located at a height of 
z = 1710 km in the middle chromosphere. The origin of 
the z axis is defined by the level where the average Rosse- 
land optical depth reaches unity. Since the horizontal ex- 
tension is 5600 km, the model represents a quadratic patch 
of 7"7 x 7."7 of a non-magnetic internetwork region like it 
should be observable with ALMA. 

The model chromosphere is dominated by shock waves 
which are excited in the lower layers. Propagation and 



interaction of these waves give the model chromosphere its 
characteristic appearance: thin filaments of hot gas (shock 
waves) and embedded cool regions (see Fig.^i). The mesh 
size of this network-like pattern is comparable to spatial 
scales typical for granulation. The average chromospheric 
temperature is roughly 3700 K but spans a large range 
from only ~ 2000 K to over 7000 K in shock waves. Caused 
by the large velocity of these waves, the whole pattern 
changes on time scales of only 20-25 s (1/e autocorrelation 
time scale). 

As mentioned in Paper I the radiative transfer in the 
simulation still needs to be improved for chromospheric 
conditions whereas the lower layers are modelled very re- 
alistically (see, e.g., Leenaarts & Wedemeyer-Bohm, sub- 
mitted). The absolute temperature amplitudes are thus 
somewhat uncertain and need to be calibrated against ob- 
servations. In contrast, the general topology of the upper 
layers (up to heights where magnetic fields become non- 
negligible) is a robust feature. The shock waves, which pro- 
duce the pattern, are generated in the lower well-modelled 
layers, and their propagation and interaction does not 
sensitively depend on the thermodynamics of the chro- 
mospheric layers. Although the temperature amplitudes 
must be considered with caution, the model nevertheless 
provides insight in the involved spatial and temporal scales 
and thus can help to define constraints for future obser- 
vations. 



4. Radiative transfer calculations 

The continuum at (sub-)millimeter wavelengths originates 
from the low and middle chromosphere and is mainly due 
to thermal free-free radiation. The assumption of local 
thermodynamic equilibrium (LTE) is valid and the source 
function is Planckian. Furthermore, hv <C k^T ', allowing 
to use the Rayleigh-Jeans approximation. Consequently, 
the source function is linearly related to the local gas tem- 
perature at these wavelengths. These conditions largely 
simplify the necessary calculations. 

The continuous opacity at mm and sub-mm wavelengths 
goes roughly proportional to the electron pressure. Carlsson & Stei 
showed that the ionisation fraction of hydrogen — be- 
coming the major electron contributor at chromospheric 
heights — exhibits substantial deviations from thermody- 
namic equilibrium conditions. The ionisation fraction of 
hydrogen cannot follow the rapid changes of the thermo- 
dynamic conditions associated with the passage of shock 
waves but tends to stay on a fixed level characteristic for 
the hot, shocked gas phase. In comparison to LTE con- 
ditions this leads to a systematic increase of the opac- 
ity, particularly in the cool, largely neutral regions. On 
average we expect a shift of the emitting layers towards 
larger heights, and an increase of intensity contrasts. The 
important deviations from thermodynamic equilibrium of 
the hydrogen ionisation have not been taken into account 
in the spectrum synthesis calculations of this preliminary 
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study but should be considered in more detailed calcula- 
tions. We nevertheless think that the relative spatial and 
temporal statistics of the intensity fluctuations is still rea- 
sonably represented by our simplified model. 

The continuum images, which are presented here (see 
Fig- CJ , were computed with the 3D LTE spectrum syn- 
thesis code LINFOR3D, which was developed by Steffen 
& Ludwig (based on the Kiel code LINFOR/LINLTE). 
For convenience we adopt the following wavelengths for 
the spectrum synthesis, covering the range accessible by 
ALMA: 0.3 mm (~ 1000 GHz), 1 mm (~ 300 GHz), and 
9 mm (~ 33 GHz). All calculations are done for disk- 
center (/Lt = 1.0). 

5. Results 

The horizontal temperature cross-section in Fig. la shows 
a pattern which is characteristic for the model chromo- 
sphere. The intensity maps in Fig. ^-d look more or less 
similar but with different contributions of dark and bright 
areas as can also be seen in the histograms in Fig. |5J The 
gas temperature distribution at z — 1000 km (panel a) ex- 
hibits two separate peaks, one at high temperatures due 
to hot shock waves and one at low values, caused by the 
cool background. The model shows a similar picture for 
all chromospheric heights above ~ 800 km where shock 
waves become clearly visible (see Paper I for more details). 
The layers below are characterised by a single peak with 
a height-dependent width. The differences in the intensity 
histograms are therefore directly caused by the mapping 
of different height ranges, from where the continuum in- 
tensity emerges, thus sampling different sets of gas tem- 
peratures. 

The height ranges can be seen in Fig. [3] in the form 
of intensity contribution functions. At 0.3 mm the maxi- 
mum intensity contribution originates from the upper pho- 
tosphere (z = 480 km), whereas the chromosphere con- 
tributes only little compared to the intensity at 1 mm. At 
that wavelength the middle chromosphere produces a sig- 
nificant fraction, while the maximum is still located in the 
very low chromosphere (z = 580 km). Consequently, the 
intensity distribution at 1 mm is much more similar to the 
gas temperature at z = 1000 km compared to A = 0.3 mm 
which mostly samples the high photosphere. 

At 9 mm, however, the intensity is due to a large height 
range throughout the whole model chromosphere, even ex- 
ceeding its vertical extent. Hence, a significant part of the 
emergent radiation at A = 9 mm would come from layers 
not included in the model, so that the results derived for 
this wavelength should be interpreted with great care. 

The clear tendency of sampling higher layers as one 
goes to longer wavelength is expected since the opacity 
scales as the wavelength squared. The rapid increase of 
the opacity towards longer wavelength is associated with a 
shift of the continuum forming layers towards correspond- 
ingly lower density. 



In principle, the temporal evolution of the chromo- 
spheric structure depends on the sampled height range, 
too. But since the evolution time scale (see Paper I) stays 
roughly constant throughout the model chromosphere, the 
intensity maps evolve on the same time scales (20-25 s). 

Finally, we degraded the spatial resolution with a Gauss 
profile. The small-scale pattern is clearly visible at 0."5 
(FWHM). At 1", however, the image is very blurred but 
still exhibits fluctuations on granular scales. But the in- 
tensity contrast decreases strongly with resolution, com- 
plicating the detection of the pattern at resolutions above 
1". Pushing the effective resolution as high as possible is 
thus clearly of crucial importance. These matters are dis- 
cussed in more detail in a forthcoming paper. 

6. Conclusions 

A comparison with future ALMA observations will form a 
strong test for present and future models of the solar atmo- 
sphere. As shown with these exemplary calculations, emer- 
gent continuum intensity at different wavelengths sam- 
ples different height ranges of the (model) chromosphere. 
Hence, simultaneous observations at multiple wavelengths 
will provide substanial new insight in the hitherto hotly 
debated structure of the solar chromosphere. Given this 
powerful tool, also the propagation of shock waves and 
chromospheric oscillation modes can be addressed. How- 
ever, as we did show in this work, high spatial and tem- 
poral resolution are of crucial importance for observing 
the small-scale structure of internetwork regions, whereas 
the field of view just needs to be large enough to cover a 
reasonable part of an internetwork region. 
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Figure 2. Histograms for a) gas temperature in a 
zontal slice at geometrical height z = 1000 km, and 
emergent continuum intensity at wavelengths of 0.3 
1 mm, and 9 mm, respectively. 
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Figure 1. Exemplary horizontal maps, all calculated from 
the same simulation snapshot: a) gas temperature at a 
geometrical height of z = 1000 km, and b-d) emergent 
continuum intensity at wavelengths of 0.3 mm, 1 mm, and 
9 mm, respectively. The legend next to each panel indicates 
data range and colour- coding. All images are calculated for 
disk-center. 
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Figure 3. Normalised averaged intensity contribution func- 
tions on the geometrical height scale of the simulation se- 
quence at wavelengths of a) 0.3 mm, b) 1 mm, and c) 
9 mm. The vertical lines mark the height of maximum 
contribution (dashed), of half the integrated function (dot- 
dashed), and the boundary photosphere/chromosphere 
(dotted). 



